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Abstract 



<d : 

. We study the consequences of imposing an approximate Galilean symmetry on the Effective Theory 
of Inflation, the theory of small perturbations around the inflationary background. This approach 
allows us to study the effect of operators with two derivatives on each field, which can be the leading 
■^j- ' interactions due to non-renormalization properties of the Galilean Lagrangian. In this case cubic non- 
Gaussianities are given by three independent operators, containing up to six derivatives, two with a 
(T) ■ shape close to equilateral and one peaking on flattened isosceles triangles. The four-point function 
is larger than in models with small speed of sound and potentially observable with the Planck satellite. 



1 Introduction 



In this paper we study the observational consequences of the presence of an approximate Galilean 
symmetry [T] in the theory of inflation, or more precisely in the theory describing the perturbations 
around an inflating solution. While it is always important to study all the possible symmetries that 
can constrain the dynamics of a physical system, this is particularly true for single field inflation 
since, given its simplicity, there are not many symmetries one can impose. In a recent paper [2] (see 
also [22 H|) Burrage et al. studied inflaton models endowed with the Galilean symmetry 

0-*-^ + 6 M a^ + c, (1) 

with b^ and c constant. They restricted their analysis to Galilean operators of the schematic form 
(d(f)) 2 (d 2 (j)) n , n < 3, which are the ones that give second order equations of motion of the form 
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(d 2 cj)) n+1 [1\ . In this way perturbations remain well behaved, even when these operators are impor- 
tant for the background solution, which is the most interesting regime, quite different from standard 
slow-roll. Interestingly this model can give rise to large non-Gaussianities and, even more interest- 
ingly, the cubic operators which generate them contain four derivatives, so that they are naively 
distinguishable from the ones generated in more standard models, where large non-Gaussianities are 
linked to a low speed of sound and the cubic operators have three derivatives Unfortunately, 
all the operators with four derivatives can be rewritten, as we will explain below, in terms of the ones 
with three using the linear equations of motion; this implies that this class of models does not give 
rise to a distinctive shape of non-Gaussianity. This fact motivates us to look beyond the Galilean 
operators studied in [2], in particular at the ones with two derivatives on each field. We will see 
that this is consistent - Galilean operators with less derivatives are not radiatively generated - and 
it gives new shapes for non-Gaussianity, with a maximum of six derivatives. 

At first the idea of studying inflationary models where operators of the form (d 2 (p) n are important 
seems silly. First of all there is an infinite number of operators of this form, with all the possible 
contractions of indices. Moreover when these operators become important for the unperturbed so- 
lution - which is mandatory if we want them to be relevant for the dynamics of the perturbations 
- the equations of motion are of higher order and describe also ghost excitations. In other words 
the solution we are looking seems outside the regime of validity of the effective field theory (EFT). 
However this comes from a wrong approach to the problem. What matters is the existence of a good 
EFT describing fluctuations around a quasi de-Sitter background: this is what we call inflation and 
this is what is relevant for observations. This is the approach advocated in [7] (see also [8]) and we 
will review it in this context in Section [2j The construction of the theory and its radiative stability 
are discussed in Section [3] and Appendix [A] The study of the three-point function is done in Sections 
13.21 and 14} while Section [5] is devoted to the four-point function. Conclusions are drawn in Section [6l 



2 The (Galilean) effective theory of inflation 

Inflationary observables - the n-point functions of the conserved curvature perturbation on comoving 
slices, which we call C 0) 

- are calculated within the theory describing small fluctuations around a 
quasi de Sitter background. As in any effective theory, we need the cut-off A to be much larger than 
the typical energy scale of the problem, i.e. H, which is the scale of quantum fluctuations induced 
by the cosmological background. This is all we need and what we call effective approach to inflation 

However one may be more ambitious and require the validity of the theory to be much broader 
than what is needed to reproduce cosmological observations. Although we experimentally probe 
small perturbations 

(t>(t + iT(t,x)) , ^ = -C~10" 5 , (2) 

one may wonder whether the regime of validity of the theory encompasses very different backgrounds: 
for example whether it can describe a solution that starts with zero velocity <p = or that, after 
the end of inflation, oscillates around a local minimum. In other words whether the theory makes 

1 This quantity is usually called 1Z. However, we follow the notation of |13j . which has become common in 
the literature. 
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sense also for tt ~ t. This is often implicitly assumed in the standard approach to inflation, where 
one starts from a Lagrangian, finds a classical inflating solution and then studies small perturbations 
around it. We stress that, although it is nice to have a single EFT describing quite different classical 
solutions, that can be used not only to derive the inflationary observables but also, for instance, to 
address the issue of initial conditions and reheating, this is not required and it may represent an 
unjustified restriction on the theory of inflation. For large values of tt the EFT may cease to make 
sense and the appropriate description may be in terms of completely different degrees of freedom. 

Usually the issue is not appreciated as in the simplest model of inflation, obtained expanding 
around a classical slow-roll solution of a Lagrangian for a single scalar with minimal kinetic term and 
potential, one can indeed consistently describe very different backgrounds, assuming a sufficiently 
heavy UV completion. Things are a bit subtler for inflationary theories with a small speed of sound 
for perturbations. Although they can be described in the effective approach [7] in a simple way, these 
theories can be also explicitly derived from a Lagrangian of the form P(X, 0), X = — {dip) 2 . If P is 
given, one can trust solutions quite different from the inflating one, at least as long as the kinetic 
term of perturbations remains healthy [9]. This of course is possible only at a price: while the theory 
of perturbations around the inflating solution just depends on few parameters, one must know the 
whole function P(X, (j>) from some UV input if interested in very different backgrounds. 

In the case under study, however, one is forced to give up the (unnecessary) ambition of describing 
solutions which are very far from the inflating one. We will be interested in a Lagrangian for 
perturbations of the schematic form 

M^H(8tt) 2 + M(d 2 Trf + ... (3) 

and particularly in the regime of large non-Gaussianities, when the second term, evaluated at d ~ H 
and Hit ~ I CP 5 , gives a non negligible correction to the kinetic term, say of order 10 -3 to be around 
the present experimental limits. If we now try to describe solutions with d ~ H and tt ~ t (which 
implies Hit 3> 1 as H ~ e/t in terms of the slow-roll parameter e = —H/H 2 ) the importance 
of the cubic term compared to the quadratic will increase by a factor 10 5 /e, boosted by the large 
classical occupation number of tt. In the EFT approach however, higher derivative terms must be 
small corrections and must be evaluated using the lowest order equations of motion (see for example 
|10[ If on the other hand higher derivative terms are not small, one should solve the complete 

equations of motion which would, in this case, contain higher time derivatives and need additional 
initial conditions. These new degrees of freedom are ghost-like, a clear sign that we are out of the 
regime of validity of the EFT we started with. 

What are the rules to build the effective theory of inflation? In [8], the emphasis was put 
on a geometrical approach. Here this would not be very useful, because it is difficult to construct 
Galilean operators in the geometrical languagd^]. Moreover we will see that the mixing with gravity 
is negligible so that one is solely interested in the action for tt. The straightforward approach is 
to write down all the possible tt operators which are Galilean invariant and realize nonlinearly the 
Lorentz symmetry. This last requirement is equivalent to build Lorentz invariant operators in terms 

2 For example one of the key point of the geometrical construction is that, given that there is a privileged 
slicing of the spacetime, one can project all tensors parallely or perpendicularly to it. However these projectors, 
once expressed in terms of the scalar degree of freedom, are not Galilean invariant as they contain a single 
derivative. 
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of ip = t + Tr(t,x), as tp linearly realizes the Lorentz symmetry. Once the possible operators at each 
order are written in terms of tt the size of each term can be estimated by naturalness arguments. 



3 Building up the action 



The purpose of this paper is to study the theory of perturbations about an inflating background when 
they are endowed with the Galilean symmetry tt — > tt + b^x^ + c; in particular we are interested in 
seeing whether we can generate new shapes for the three-point function. Galilean invariance implies 
that there are at least two derivatives for each tt in the equation of motion, which means that a cubic 
term in the action contains at least four derivatives (we are assuming that we can take the decoupling 
limit and just look at the action for tt, neglecting the mixing with gravity; we will come back to this 
point later). One has cubic terms with only four derivatives if the action for tt is completed to 
Lorentz invariant operators for ip that have the minimum number of derivatives compatible with the 
Galilean symmetry, that is (dil)) 2 (d 2 ip) n , because these terms give equations of motion with exactly 
two derivatives on each tt [1]. 

It is easy to see, however, that cubic operators for tt with four derivatives do not generate new 
shapes of non-Gaussianity as they can be rewritten in terms of operators with three derivatives, 
which are the ones typical of models with a reduced speed of sound^ [6j [7]. Indeed the possible 
forms of the four-derivative operators are 



where with V we indicate spatial derivative^. To show that they can be rewritten as operators 
with fewer derivatives, we are going to impose the linear equation of motion (at leading order in 
slow-roll and with generic speed of sound c s ) tt + 3Htt — c 2 X7 2 Tr/a 2 = at the operator level; this 
is equivalent up to third order to perform a field redefinition. Since we want to compute equal-time 
correlation functions and not S-matrix elements, in principle one should take into account all the field 
redefinitions [13] . However it is easy to realize that the redefinition we need involves derivatives of the 
field and therefore, when the modes are much longer than the Hubble radius, it gives contributions 
of higher order in the slow-roll parameters. 

The first operator in @ is obviously a total derivative, so that a time derivative can be moved 
to act on the a 3 of y/— g giving a term with three derivatives. The second is related to the first using 
the linear equation of motion, apart from terms with fewer derivatives. The third one can be written 
as ^VttVtt 2 and it is thus related to the second integrating by parts. The fourth is the same as the 
third integrating by parts in time. And finally the fifth is related to the fourth using again the linear 
equation of motion. Therefore with four derivatives we cannot generate new forms of three-point 
function and we are hence motivated to consider terms with more derivatives on ip. 

3 We thank A. Gruzinov, M. Mirbabayi, L. Senatore and M. Zaldarriaga for pointing this out to us. 

4 Reference [3] claims that the shape induced by the fifth operator in (TjTJ is independent of the shapes given 
by operators with three derivatives. However it is straightforward to check that it can be written as a linear 
combination of the shapes induced by tt 3 and 7r(V7r) 2 . 




(4) 
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3.1 Galilean action with more derivatives 



Let us assume to start, besides the standard kinetic term, only with operators that have at least two 
derivatives on each tt (and therefore on ip), i.e. setting to zero the lowest derivative Galilean terms 
(d r ip) 2 (d 2 ip) n In flat space this is consistent as loops will not generate them |12| . However, the 
Galilean symmetry cannot be defined in curved space, since there is no natural way of defining a 
constant vector b^. Thus one may worry to generate terms containing the Riemann tensor, of the 
schematic form (dijj) 2 (d 2 il)) n R and in general all terms in which a couple of derivatives is replaced 
by the curvature tensor. All these would give a contribution to the tt correlation functions which is 
of the same order as the one we started with. This is not the case as it is easy to prove. 

Let us start from a 1PI graph for tt in flat space: this has (at least) two derivatives on each 
external leg. Now we add an external graviton linfi If it couples to an internal line, it will not 
reduce the number of derivatives acting on each external tt line. If the graviton is attached to a 
vertex it can reduce to one the number of derivatives acting on an external tt line, when g comes 
from the Christoffel symbol. But all the other external tt legs have still 2 derivatives. The loop can 
generate additional external derivatives, but it cannot reduce them, so we see that it is impossible 
to generate an operator of the form (dip) 2 (d 2 ip) n R, which contains a term with a single external 
graviton and two 7r's with a single derivative. One would have at least an additional suppression of 
the curvature: (dip) 2 (d 2 ip) n R 2 , which makes these terms subleading for the tt correlation functions. 
Additional non-renormalization properties are discussed in Appendix [Al We conclude that, although 
the Galilean symmetry cannot be defined in curved space, we can consistently reduce to operators 
which contain two derivatives on each field. 

Let us come to the explicit construction of the tt action. The tt operators at any order can be 
written starting from Lorentz invariant operators for tp. These are traces of the double derivative 
matrix V In the following we use ^ to denote this matrix and the square brackets to indicate 
the trace so that, for example, [tyty] means V^V^ V^V 1 '^- I n order to parametrize the theory of 
perturbations around a given solution we would like to isolate operators that contain terms linear in 
tt. Those operators change the background and therefore their coefficients are constrained once the 
background solution is chosen. On the other hand operators that start quadratic in tt do not affect 
the background solution for if; nor for the metric - as they have vanishing stress-energy tensor for 
tt = - so that their coefficients are not constrained [8j [7] . 

If a term is composed by the product of several traces of products of ty, [^^f . . . \P] . . . . . . 

then we can subtract to each trace its background value to make the term at least quadratic. If on 
the contrary an operator is composed by a single trace [& n ] this cannot be done and the operator 
will start with a linear term. We will show however that [^ n ] can be written in term of operators, 
whose linear term is easy to isolate. To do this, consider the sum 

£(_ . . . g^^V^V^ . . . V^V,^ , (5) 
p 

where p is the parity of the permutation. This sum contains also the single trace of order n. For n > 4 
eq. ([5|) trivially vanishes, as there are too many indices to antisymmetrize, so that the [^ ra ] can be 
rewritten in terms of products of two or more shorter traces. For n < 4 the sum does not vanish. In 

5 Notice that this implies that the kinetic term is standard, i.e. c s = 1. 

6 We do not consider internal graviton lines, as these would be suppressed by powers of Mpi. 
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flat space it is a total derivative, as it easy to see calculating the equations of motion, which trivially 
vanish, but in curved space this does not occur since derivatives do not commute. Nevertheless the 
sum can be written — integrating by parts one of its derivatives, writing the commutator of derivatives 
in terms of the Riemann tensor and finally expressing the Riemann tensor of de Sitter in terms of 
the metric — as 

H 2 ^(-l)yi^) . . . ^-xpK-Ov^V^ • • • V^-iV^-i^ n < 4 , (6) 
p 

apart from overall combinatorial factors. Notice that these are the minimal Galilean terms, see 
Appendix A of [T], up to the quartic. Therefore for n < 4 the single traces can be written in terms 
of shorter traces and minimal Galilean terms, cubic and quartic in tp (0). We conclude that it is 
consistent to study the theory that contains all operators written in terms of V^V^i except the 
single trace operators, plus the minimal Galilean operators of cubic and quartic order. These two 
will be suppressed by H 2 and therefore will give a contribution of the same order of magnitude as 
the other ones to the 7r correlation functions. 

It is straightforward to write the minimal Galilean operators of cubic and quartic order in such a 
way that they do not contain terms linear in ir. The cubic (DGP-like) term Otp(dip) 2 can be written 
as (Dtp + 3H) [(dtp) 2 + 1], just redefining the coefficients of the kinetic term and the cosmological 
constant term, which are anyway present in the Lagrangian. The procedure is slightly subtler for 
the quartic term and it is useful to rewrite it as |14| 

(aip) 2 (d%b) 2 - 2nv;d At ^v At v^ - (v M v^) 2 (9V) 2 + 2d^v^v^v u v p ^d p ?p 

=(«9V0 2 [2(D^) 2 - 2(V a V / #) 2 ] - (d^R^d^d^ . 
As we are in de Sitter the Ricci tensor is proportional to the metric, so that the second operator is 
just (dip) 2 (dtp) 2 . In this form all terms can be treated similarly to what we did with the DGP one. 

Now that we have the rules to build the action in terms of ip, we move to study the third order 
action for ir, which is relevant for the calculation of non-Gaussianity. 



3.2 Cubic action for perturbations 

We are primarily interested in the 3-point function, i.e. in operators of the form (d 2 Tr) 3 . Some of these 
operators will be independent, while others will be related to quadratic operators by the non-linear 
realization of the Lorentz symmetry, i.e. once written in terms of ip. In this second case the scale 
suppressing the cubic operators is related to the one suppressing the quadratic ones. Schematically 
we will have 

M^H(dir) 2 + M [H(d\) 2 + (d 2 irf + . . .] , (8) 

where the factor of H inside brackets comes from the fact that V^V^V gives H5ij on the background. 
Going to canonical normalization we have 

(dTr c ) 2 + ± I (d 2 Tr c ) 2 + j- 5 (d 2 ir c f + ... (9) 

7 This represents a nice consistency check of our statement before, that generic operators of the form 
(dip) 2 (d 2 ip) n R are not generated. As the operator ^ is Galilean invariant, it will not generate non-Galilean 
terms, unless paying an additional power of R, i.e. (dip) 2 (d 2 ip) n R 2 , consistently with what we said above. See 
also Appendix [21 Moreover, as minimal Galilean terms are not renormalized, we are not going to generate 
the quintic one, which does not come in the process of rewriting the single trace operators. 
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where 

.2 M ^H , . 

~ ~HM~ ( 10 > 

and 

A = A — ~ A — • \fe — — fill 

1V J1 ^3 A 3 H 

where the separation between the scales is stable only when A > A. Notice that (<9 2 7r) 2 must be 
treated as a small perturbation to the standard kinetic term with c s = 1; therefore we have the 
usual normalization H / (-y/eMpi) ~ 10~ 5 . Cubic non-Gaussianities will be of order (H/A) 5 , so that 
observable non-Gaussianity - i.e. larger than 10 -5 corresponding, in the usual parametrization, to 
/nl ^ 1 ~~ implies A> A. We conclude that, assuming all independent operators are suppressed by 
a common scale A, cubic operators that are linked by symmetry to quadratic ones are suppressed by 
a much higher scale and can be neglected with respect to the independent ones. It is straightforward 
to see that this holds in general when an operator (d 2 ir) n is related by the Lorentz symmetry to a 
lower dimensional one: the one with the lowest dimension will be suppressed by the common scale 
A, while the others will be suppressed by powers of ^ ■ -y/e^p. The same conclusion applies to the 
minimal Galilean operators of cubic and quartic order, where there is a relation between quadratic 
and cubic terms in ir. Therefore these operators have a negligible contribution to the three-point 
functions in our setup. 

Since we are interested in quadratic and cubic operators and we do not have linear terms, we 
need to consider only terms with two traces on tjj and terms with three traces: 

([**... ¥]-ci)([¥¥...¥]-c&) and ([**... ... *]-c 4 )([** ... *]-<%) (12) 

where the constants Cj are chosen to subtract from each trace its background value. Operators that 
contain a Uir are proportional to the linear equation of motion and can be neglected since they can 
be removed by a field redefinition, as we discussed above. Of course this is valid only at leading order 
in slow-roll, but still tells us that operators with a Ott are subdominant with respect to the others. 
The only cubic operator without Uir, coming from a term with three traces is of course 

(V 2 tt) 3 V 2 = <f J 'ViVjTr . (13) 

Terms with two traces generate terms with 2 and 3 7r's and it is important to understand whether the 
latter are independent or not because this controls the scale that suppresses them, as we discussed 
in the previous paragraphs. If we restrict again to cubic operators without Utt the additional terms 
are 

V 2 7r(V,V^) 2 V 2 vr(V i V / ,^) 2 vMv^tt) 2 . (14) 

It is straightforward to see that there is enough freedom to make all of them independent from the 
2 7r's operator^. 

So far we have implicitly assumed that we can concentrate on the ir action, without taking into 
account its mixing with gravity. To check that this mixing is indeed negligible, one should follow 



8 This independence is not general: not all the cubic operators (9 2 7r) 3 are independent from the quadratic 
ones. It is straightforward, however, to check that the ones which are dependent contain a Ott and are thus 
neglected here. Another example is given by the two minimal Galilean operators that we are considering: also 
these are thus neglected in the calculation of the three-point function. 
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in spatially flat gauge one solves the ADM constraints and plug the solution for the ADM variable ./V 
and Ni back into the action. To derive the cubic action, one needs to solve for ./V and 2Vj at first order 
only, as the second order result would multiply in the action the first order constraint equations, which 
vanish |13| . If the quadratic action is simply given by the standard kinetic term, then the mixing 
with gravity at energies of order H is suppressed by the slow-roll parameters, so that the mixing with 
gravity just gives corrections to the 3-point function of the order of the slow-roll parameters. In our 
case we also have additional quadratic operators - see eq. ([9]) - of the schematic form (d 2 TT c ) 2 / A 2 . 
However, neglecting the mixing with gravity, these terms are suppressed with respect to the canonical 
kinetic term by H 2 /A 2 , so that we expect they also give a suppressed contribution to the constraint 
equations. Indeed it is straigthforward to check that their effect is suppressed both by slow- roll and 
by H 2 /A 2 . Therefore we can safely concentrate on the ir action, without bothering about the mixing 
with gravity, to calculate the observable predictions. 

At third order in perturbations, we have the following four independent operators: 



S 



M^Hd^TT + Mi^ViVj-vr) 3 + M 2 (5 ii V i V i 7r)(V i V i 7r) 



+ M z (g ij ViV^){ViV^) 2 + M 4 (5 4j V i V i 7r)(V / ,V,7r) 2 



in which we need to expand the covariant derivatives: 



fl ij ViV,-7T 



o) 



7T 



a 2 u V 



d 2 TT 



3Hn, 



(15) 



(16) 



{ViVjiry = — j<5 <5 ((?jSj7r - Ha^^idkdnr - Ha z 5 kl k) 



(V^Vovr) 2 = - — (diTT-Hdiir) 



2Htt-V + , (17) 

a z 

(18) 



As discussed in the previous subsection, in order to calculate the non-Gaussian correlation functions 
at lowest order in slow-roll one can substitute the linear equation of motion jt + 3Hn = V 2 7r/a 2 in 
the action. Substituting the Laplacian in favour of the time derivatives (in order to recover the usual 
shapes involving the operators fr 3 and 7r(<9j7r) 2 /a 2 ) yields 



gVViVjTT = 7T , (ViVj 



7T 



2Hnjf - m 2 -k 2 



(19) 



which results in an action with only three, and not four, independent operators. Doing a further 
integration by parts, a^Tck 2 = ^a 3 dtTT 3 — > —a 3 H7r 3 , we finally get 



S 



dW 



M^H ( -k 2 - ^jt- ) + Miff 3 + M 2 ff- 



+ M 3 7f 



.. (didjTr) 



2H-frjt 2 + 3i? 3 vr 3 



(20) 



This will be the starting point for the calculation of the three-point function. 



4 Shapes of non-Gaussianity 

Using the cubic action we can calculate the tree-level contribution to the three-point correlation 
function for the tt field. With the prescriptions of the in-in formalism, this is given by the following 
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formula [T3j [15] : 



(- 4 (/u)) - / r dta 3 ( 



7T 3 (i ), / d 3 xL int (i,x) 



(21) 



where to is some late time when the physical modes of interest are well outside the Hubble radius, 
and Li n t is the cubic Lagrangian to be read from Eq. (|20p . 

Non-Gaussianity is usually computed in terms of the correlation functions of the Bardeen potential 
<£. It is customary to write the correlation functions for $ isolating the Dirac's delta that enforces 
the conservation of the three-momentum (implied by the traslational invariance). In the limit of 
exact scale invariance, the two- and three-point functions can be parametrized as 

<%A> = (2^(k 1 + k 2 )^ (22) 

and 

= (2vr) 3 ,5( 3 )(fc 1 + k 2 + k 3 )F(h, k 2 ,k 3 ) , (23) 

where F (called the bispectrum of <£) is a homogeneous function of degree -6 which, because of 
rotational invariance, depends only on the length of the three-momenta. Therefore, in the same limit 
kfF(k\,k 2 ,k 3 ) becomes just a function of the ratios r 2 = k 2 jk\ and r 3 = k 3 jk\. As argued in |16j . 
the most meaningful quantity to plot is 

rir 2 3 F{l,r 2 ,r 3 ) , (24) 

which is also commonly called the "shape" of the non-Gaussianity. 

During matter dominance <3? is related to the curvature perturbation on comoving slices £ by the 
simple relation $ = |(\ The curvature perturbation is constant outside the Hubble radius, regardless 
of the matter content of the Universe (and in particular it remains constant during reheating). At 
leading order in slow roll £ is related to the tt field by £ = —Hit. Shifting to conformal time 
(r = —1/aH), the bispectrum F reads thus 

i /3\ 3 f° dr 



F(k 1: k 2 ,k 3 ) = -—^-j 7r*(0,fci)7r*(0,fc2)7r*(0,fc3) J L int (r, k u k 2 , k 3 ) + c.c. + perms. , (25) 



where 7t(t, k) is the classical Fourier mode in de Sitter space (see Appendix [B]) and the momenta are 
restricted to configurations forming a triangle in momentum space. 

Doing this for each of the three independent cubic interactions in (|20p leads to three different 
bispectra, which are proportional to Mi, M 2 and M3 respectively. As detailed in App. [Bl their form 
is (setting kt = k\ + k 2 + k 3 ) 

„ fu , ,x 20 M X H 2 kf - 6kt(hk 2 + k 2 k 3 + hk 3 ) + 15hk 2 k 3 

F Ml ( kl ,k 2 , h)=~ Y7 ^ 6A $ , (26) 

5 M 2 H 2 k\ - k\ - k\ \ t 2{k 2 + k 3 ) - h 



F M2 (h,k 2 ,h) =6^A,^ pp< 



2 + 



n h{k 2 + fc 3 ) + 2k\ - 2k 2 k 3 + 2fc| h(j4 - k 2 k 3 + k\) - 2k\k 3 - 2k 2 k 2 3 
k 2 + k 3 



2i kik 2 k 3 {k 2 + k 3 ) + 2klkl _ fcifc|fc§ 



+ cyclic permutations , (27) 
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FM 3 (k 1 ,k 2 ,k 3 



5 M 3 H 

3 eMl kik 2 k 3 k? 



A 2 



+ 



'PI 

k\ 



ko + A;? ^3 

10-24 J + 48- 



&2&3 

+ cyclic permutations . 



1 + 3 



&2 + k 3 



h/2 k x {k 2 + k 3 ) - 2k 2 k 3 
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30 



^^2^3 



(28) 



The simplest shape of non-Gaussianity is the so-called local shape, which peaks in the squeezed 
limit (r 2 — 1 and r 3 ~ 0, i. e. when a mode is much larger than the other two). In single field inflation 
non-Gaussianity of this kind is proportional to the tilt of the power spectrum, and a measurable signal 
can be generated only within multi- field scenarios |23j. Single field inflation can generate large non- 
Gaussianity when derivative interactions are present, whose shape is enhanced in other regions of the 
r i — r 2 plane than the local limit. The shape generated in scenarios like DBI inflations, for instance, 
peaks in the equilateral limit (r 2 — ^3 ~ 1) [16J. Shapes peaking in the "flat" configuration (r 2 ~ 
r 3 ~ 0.5, corresponding to flattened isosceles triangles) where obtained in Ref. |17| as a difference 
of equilateral shapes, and directly from new higher derivative operators in Ref. |18| . Notice however 
that Ref. [18J did not identify any symmetry that allows to neglect lower dimensional operators and 
that operators containing higher time derivatives were not considered. 

The shapes given by the three independent cubic operators in Eq. (|20p are shown in Fig. [TJ The 
first operator 7f 3 has a clearly equilateral shape and is almost indistinguishable from the operator 7r 3 : 
indeed, in conformal time 7^ ~ k 2 T 2 e tkr is very similar to tt^ ~ k 2 T 2 {l + ikr /2)e tkT . For the same 
reason, 7f(<9j7T — Hdiit) 2 is very close 7r((9j7r) 2 . Both operators have therefore equilateral shapes of the 
kind already considered in Refs. [5j[6lE]- 

On the other hand, an asymmetric operator with more scalar product of space gradients like 
7f(9j3j7r) 2 /a 4 will contain terms like (k 2 ■ k 3 ) 2 = \{k\ — k\ — fe|) 2 - Moreover, each power of 1/a ~ r 
yields 1/(^1 + ^2 + ^3) after the integral over conformal time; scalar products of gradients will therefore 
result in terms like [(k 2 — k 2 — k 2 )/2{k\ + k 2 + k 3 ) 2 ] n , which tend to be suppressed in the equilateral 
limit. The precise scaling depends on the details of the symmetrization, but the general argument 
qualitatively holds. As a result, the shape generated by the operator proportional to M 3 looks like 
a wave reaching its maximum in the flat configuration. 

In order to give a quantitative measure of the difference of the shapes generated by different 
interactions, it is very useful to define a scalar product between bispectra as follows |16] : 

„ „ F(i){ki,k 2 ,h)F {2) (ki,k 2 ,k 3 ) 

" ^ = f A, P(h)P(h2)P(h) ^ 

triangles 

where the sum is restricted to the vectors hi that form a triangle in momentum spac^§, from which 
it is natural to define the cosine 

cos(F (1) ,F (2) ) = - F F § 2 ^ 2) F )l/2 ■ (30) 

If the cosine between two shapes is large, they have a significant overlap and it is possible to use 
the same estimator to analyze CMB data and constrain the amplitude /nl of each shape. On the 



For CMB applications it is useful to introduce a 2d scalar product that quantifies the similarity of two 
shapes in a given CMB experiment |TB] . Here we stick to the 3d definition. 
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Figure 1: Bispectrum shapes generated by the cubic operators in Eq. (|20p . The ones proportional to 
Mi and M2 (upper panels) are very similar (with a cosine of nearly 1) to those generated by 7r 3 and 
7r(<9j7r) 2 respectively, which where studied in [5j El [7]. They have a cosine of 0.96 and 0.99 with the 
equilateral template. The third operator (lower left panel) is "surfing", i.e. it is shaped like a wave 
with a maximum in the flat configuration, and has a large overlap with the enfolded template (lower 
right panel). The suppression in the equilateral regime is due to the presence of a higher number of 
scalar products of gradients. 

other hand, if the cosine is small this is no longer efficient and a different estimator is needed in order 
to set optimal constraints on /nl- 

Moreover, a crucial numerical boost in CMB data analysis is gained by using factorizable shapes, 
i.e. shapes which can be written as sums of monomials of k\, k% and ^3 |19|, 120], A good estimator 
to constrain the amplitude of a given model is then provided by a factorizable shape having a 
large cosine with the bispectrum generated by the model itself. Factorizable shapes which resemble 
model predictions are often called templates. Among the most popular templates are the local 
and equilateral ones |16j . for which CMB constraints are usually cited (see for example |21j). Linear 
combinations of equilateral operators can however yield a significantly different shape for a quite wide 
range of coefficients, leading to the definition of an "orthogonal" template |17| . Another template 
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Shape: 


local 


equilateral 


orthogonal 


enfolded 


Ml 


0.49 


0.96 


-0.06 


0.70 


M2 


-0.43 


-0.99 


-0.16 


-0.53 


M3 


0.57 


0.71 


-0.52 


0.94 



Table 1: Cosines between different shapes of bispectra. 



found in the literature is the one for the so-called "enfolded" shape, originally introduced in Ref. |22| . 
This template is a linear combination of the orthogonal and equilateral templates, as pointed out in 
Ref. [IT]: 



1 



F cni (h,k 2 , k 3 ) = - F c ^ n (h,k 2 , k 3 ) - F mth (k u k 2 , k 3 ) 



1 1 

+ 



1,3 J>3 h,2h,2h,2 hilc 2 !? 3 lcr,k 2 k 3 
12 12 3 1 2 3 tv^rv^riiQ 



+ cyclic permutations 



(31) 



We have computed the cosines of our shapes with the local, equilateral, orthogonal and enfolded 
templates. Our results are summarized in Table [JJ As seen in this table, a suitable template for the 
operators 7r 3 and 7f(<9j7r — Hd{Tr) 2 is the equilateral one. The remaining operator has a large overlap 
with the enfolded template@. Comparing with the templates allows us to compute the amplitude of 
non-Gaussianity of each type generated within our model, and hence constrain Mi, M 2 and M 3 . 

For nearly equilateral shapes, the amplitude of the non-Gaussianity is defined as _/S = k e F(k, k, k) /6A| 
[19J. The amplitude induced by the two shapes Fm\ and Fm 2 is thus 



! = _80_ M X H 
/nl 729 eMlj 



using the most up-to date values /, eq 



NL 



26 ± 140 



240 ± 1280 



,eq,2 = __865_ M 2 H 
/nl ' 2916 £M|j 



CL) from Ref. [21 J one gets 
M 2 H 



(32) 



eMl x 



-80 ± 470 



(33) 



For the shape Fm 3 the amplitude must be defined starting from the enfolded template, 

k 6 F Ma (k,k/2,k/2) 35 M 3 H 



/•enf ,3 
•/NL 



96A| 



256 eM|j 



(34) 



Since this template is a linear combination of the equilateral and orthogonal ones (which are slightly 



correlated), the errors on the value of /§l 

can be inferred from the errors on ^^d given 
in Ref. |21j . assuming that the correlation coefficient remained the same as in Ref. [1 
obtained in this way is 

/§ n L f = 114 ± 72 (68% CL) , 



The limit 



(35) 



10 There exists a linear combination of the equilateral and orthogonal templates, F oc — 0.5i ;iorth + 0.86i ;icqul1 , 
which has an even larger cosine (0.98) with Fm 3 - However, it is good enough for our purpose to stick to the 
enfolded template, which has a cosine of 0.94 and is already known in the literature. 
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leading to 

M 3 H , . 

—V = 830 ± 530 . 36) 

A drawback of the enfolded template (as well as of the orthogonal one) is that it tends to a 
constant in the squeezed limit, while in the approximation of exact scale invariance any single- 
field three-point function vanishes in this regime [23J. This means that, while acceptable for CMB 
data analysis, such templates cannot give satisfactory results in contexts where all the signal is 
concentrated in the squeezed configurations as, for instance, with halo bias observations |24 [ I25 [ [26] . 
In order to improve the above constraint by combining it with LSS observations is therefore important 
to devise a different template, that peaks in the flattened limit but preserves the correct behavior 
when k\ <C ~ k^. An example of such a template is given in App. O 

In our theory there are 3 independent parameters which allow us to choose arbitrary linear 
combinations of our operators. Fixing the amplitude of non-Gaussianity leaves us with a 2-parameter 
family [MijM\ and M3/M1) of different possible shapes. It is interesting to see whether one can 
generate a suitable linear combination that is nearly orthogonal to all of the existing templates, 
following the same reasoning that lead to the definition of the orthogonal template in Ref. |17| . If 
this were the case, one should use a dedicated template for data analysis. 

If such a linear combination exists, the vector (1, M2/M1, M3/M1) must be an approximate null 
eigenvector of the 3x3 matrix of the scalar products of the three shapes (with the Mi's set to 1) with 
the local, equilateral and orthogonal templates (the enfolded is already a linear combination of two 
of them). Equivalently, the coefficients (Mi/Mi)\\FMi || must form an approximate null eigenvector 
of the 3x3 matrix of the cosines, which are just normalized scalar products. Since the determinant 
of the matrix is 0.0057 (much smaller than most of the entries), such approximate eigenvector can 
be found and this allows to solve for M2/M1 and M3/M1. For instance, a linear combination with 
M2 = 0.32Mi and M3 = — 0.42Mi is a good candidate, yielding cosines of -0.15, 0.03 and 0.06 with 
the local, equilateral and orthogonal template (and -0.03 with the enfolded template). The values of 
M2/M1 and M3/M1 yielding small (i.e. < 0.2) cosines with each template are shown in Fig. [2l where 
the shape of the combination nearly orthogonal to all templates is also given. 



5 Four-point function 

Let us see what are the implications of our setup for the four-point functional. The action, besides 
the quadratic and cubic operators we considered, will contain quartic operators 

(^c) 2 + ^(d\ c f + ^{3\ c f + ^(dV c ) 4 + . . . (37) 

The deviation from Gaussianity induced by the quartic terms can be quantified comparing these with 
the kinetic term at scales of order H 



G4 = W Pi 



8 



f ] 



The contribution to the four-point function of operators with higher spatial derivatives was considered in 
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Figure 2: Left: regions of the parameter space M2/M1 and M3/M1 for which the cosine of the 
resulting shape with the local, equilateral and orthogonal template is smaller than 0.2. Right: Non- 
Gaussian shape obtained with Mi = 0.32Mi and M3 = — 0.42Mi. This shape is nearly orthogonal to 
all known templates (with cosines of -0.15, 0.03, 0.06 and -0.03 with the local, equilateral, orthogonal 
and enfolded template), and would require a dedicated template for data analysis. 



The quartic non-Gaussianity comes from operators of dimension 12, and it is therefore suppressed 
by (if/A) 8 . Notice that in terms of the parameter tnl used in the literature, we have NG4 ~ tnlA^, 
where ~ 2 x 10~ 9 gives the normalization of the scalar power spectrum. As we saw, for the 
three-point function we have analogously 



NG 3 = ^~^ 

3 - (c 2 ) 3 / 2 ~ c 2 



X ' : m 



1/2 

indeed cubic operators are of dimension 9 and NG3 ~ /nlA^ . Without additional ingredients, we 



see from eq.s (|38j) and (|39|) that the quartic non-Gaussianity is always smaller than the cubic one 

NG 4 ~ NGg /5 , (40) 

so that the three-point function is much easier to detect. Notice, however, that the scaling of the 
quartic and cubic non-Gaussianities is different from the one one has in models with a low speed of 
sound, where the three-point function is given by operators of dimension 6 (schematically (<97r c ) 3 /A 2 ), 
while the four-point function is given by operators of dimension 8 (of the form (chr c ) 4 /A 4 ), so that 

NG3 ~ (x) NG4 ~ (x) ^ N ° 4 ~ N ° 3 ' ^ 

This implies that, at the same level of cubic non-Gaussianity, our setup gives a larger four-point 
function signal, although still much smaller than the cubic one. This difference is somewhat important 
experimentally. If the three-point function is of the order of the present experimental limit, /nl ~ 
100, equivalent to NG3 ~ 5x 10 -3 , then in models with a small speed of sound one has NG4 ~ 2x 10 -5 , 
i.e. tnl ~ 10 4 , so that the four-point function is not detectable by the Planck satellite, while it may 
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be detectable in our setup, where NG4 ~ 2 x 10~ 4 , i.e. tnl ~ 10 5 . Obviously these are just rough 
estimates as a proper forecast for the Planck capability on the 4-point function is still lacking. 

The three-point function is usually considered the leading signal of non-Gaussianity as it typically 
dominates the higher order correlations as we saw above. However, it has been pointed out in [28j that 
it is possible to have technically natural theories in which the leading source of non-Gaussianity is the 
four-point function, as a consequence of an approximate ir — > —tt symmetry. We can straightforwardly 
apply the same arguments in our context. Let us consider the case in which, besides the standard 
kinetic term, we only have operators that start quartic in tt. Of course both the kinetic term and these 
quartic interactions respect the symmetry tt — > — tt. However, as a consequence of the non-linear 
realization of the Lorentz symmetry, these operators will also contain terms of higher order in tt, in 
particular quintic, that do not respect the synnnetrfl Obviously if these terms were unsuppressed 
the symmetry would be of no use. However, as we discussed above, operators with a given dimension 
that are related to lower dimensional ones by the Lorentz symmetry are not weighed by the cut-off 
scale A, but they are further suppressed. In particular the quintic operators will be of the form 

t<™ f ,A^.^«!. (42) 

The breaking of the tt — > —tt symmetry will thus be weighted by the small parameter ^. In particular 
loops will induce cubic operators of the form 

These radiatively induced cubic interactions only give a completely negligible 

NG -<f) 5 -(f) ! M^=(f) 2A " 2 - <44) 

We conclude that it is possible in our setup to have a large four-point function, of order (H/A) s , 
with a negligible three-point function signal. The number of independent operators - the quartic 
operators that are independent of the cubic ones - is quite large, even if we restrict to terms without 
On. Their study is beyond the scope of this paper. 



6 Conclusions and future directions 

In this paper we explored a new class of single field inflationary models, endowed with an approximate 
Galilean symmetry. Their phenomenology is quite distinctive, with new shapes of the three-point 
function and potentially large four-point function. 

The study of this class of theories can be extended in many directions. As we discussed, there 
are regions of parameter space in which the shape of the three-point function is quite different from 
all the ones for which a data analysis has been performed, so that the constraints on these regions 
are very loose. It would be important to extend the analysis of CMB maps to these shapes. Even 

12 Actually the operator (□?/; + 3H) 4 does not contain quintic or higher order terms, so that the parity 
symmetry would be exact. However this operator gives a four-point function which is four times suppressed 
by slow-roll as all the legs are proportional to the linear equation of motion. 
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more general shapes can be obtained considering, in addition to the (<9 2 7r) 3 operators, operators with 
less derivatives (Galilean invariant or not). It is technically natural to keep these additional terms 
small, since they are not radiatively induced by the operators already included in the Lagrangian. 
Therefore they can give a contribution to the three-point function comparable to that of terms with 
more derivatives. In this paper we have just performed a quite preliminary and qualitative study 
of the four-point function signal. Since in these models this signal is larger than in models with a 
reduced speed of sound (for similar values of the three-point function) and it can become very large 
imposing an approximate Z<i symmetry, we believe it is worthwhile to characterize the four-point 
signal and compare it with data. 

From the point of view of model building, it is interesting to note that one may go a step further 
and study a theory in which interactions start with three or more derivatives per field. Although this 
cannot be dictated by any symmetry - the standard kinetic term is not invariant under a constant 
shift of the second derivative - terms with less derivatives in the interactions will not be generated. It 
may be interesting to study the phenomenology of these higher derivative models. Another possible 
direction of further study is the multi-field scenario: as in the case of models with reduced speed of 
sound |29| . multi- field models may have interesting signatures and it may be useful to explore the 
effects of the Galilean symmetry in the effective theory of multi-field inflation |30j. 

Ackowledgements 

It is a pleasure to thank Andrei Gruzinov, Lam Hui, Mehrdad Mirbabayi, Alberto Nicolis and 
Matias Zaldarriaga for useful discussions. P. C. thanks Leonardo Senatore and Matias Zaldarriaga 
for stressing the effective field theory approach to inflation. The work of G. D'A. is supported by a 
James Arthur fellowship. 

A Non-renormalization of Galilean operators 

In this Appendix we complete the analysis of the non-renormalization properties of Galilean interac- 
tions in the presence of gravity. In Section [3J we discussed what happens when the leading Galilean 
terms of the form (d^)) 2 (d 2 'i\)) n are zero. In flat space these operators are not renormalized, but 
one may wonder whether in the presence of an external gravitational field generic terms of the form 
((9?/>) 4 ((9 2 ?/>) n_2 i? can be generated. This would be relevant in our construction, since Galilean terms 
with the minimum number of derivatives were introduced to get rid of terms linear in ir. We show in 
the following that operators of the form (dil)) A {d 2 ip) n ~ 2 R cannot appear, because vertices generated 
by loop diagrams with n fields if) must have at least In derivatives also on curved backgrounds. 

To make clear the logic of our argument let us review first why this is the case in flat space. We 
consider 1PI graphs that involve the leading Galilean interactions. Each vertex has at least 2 internal 
lines, otherwise the graph is not 1PI. We want to show that we can always take the two lines with 
only one derivative to be the internal ones, so that the resulting graph has more derivatives than the 
minimal Galilean vertex. Following [31], if we pick two internal lines and call them (f> the vertices 
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can be schematically of the form 

{d 2 i>) n ~ 2 d(t)d(t) (45) 

d^{d 2 ^) n ^ d 2 j>d<p (46) 

d^dil){d 2 ^) n - i d 2 <t)d 2 <t> (47) 

however it is easy to see that we can always rewrite the last two in the form of the first, namely with 
only one derivative acting on each (f). For (|46p it is immediate since 

dfj,d v (/)d a (j) = d^{d v -)4>d a 4>) - -d a (d^d v ^)) (48) 

and integrating by parts we can move one derivative to the other lines. The interaction (|47p can 
be also be put in the form (|46p . If we write it as A^ vpo d^d v (j) d p d a (j), we see that it gives a contri- 
bution to the equation of motion of the form A pvpa d^d u d p d a (j) that must vanish because the e.o.m. 
obtained from minimal Galilean terms has 2 derivatives per field. This is possible only if the totally 
symmetric part of A pvpa is zero. However, since the structure d^d u (j) d p d a (fr has less symmetries 
under permutation of indices, it doesn't imply that the vertex equally vanishes. In particular, A 
can be antisymmetric in the exchange [fip] and symmetric in all the others, thus giving a vanishing 
contribution to the e.o.m. but a non-zero vertex. In this case anyhow we can write 

d v d [p 4> d p] d a 4> = d [p (d u (f> d p] d a <j)) (49) 

and integrating by parts we end up with a vertex of the form (|46|) and then use again (|48p to complete 
the argument. 

Now we want to see what changes if we include an external gravitational field. The starting point 
is the minimal covariantization of the leading Galilean interactions and we consider again 1PI graphs 
with n external fields ip plus external graviton lines. There are two possibilities: when gravitons are 
attached to internal lines the vertices are the same as in flat space and we can proceed as before in 
order to finish with one derivative on each of the two 4>. This will simply add g or R or Vi? etc. to 
a (d 2 ip) n graph. If instead the gravitons lines go in the vertices we can try to repeat the same steps 
we did before with the substitution of ordinary derivatives with covariant ones. The extra derivative 
in the structure VV(/>V</> can be moved to the other lines using the identity ([48|) that is valid also 
with V instead of d since it involves at most two derivatives on the same scalar field and so there is 
no problem in commuting covariant derivatives. The difference arises when we consider the vertex 
A pupa V u (fi VpV CT (/>. Because A contains at most two covariant derivatives on -0 and no R, its 
symmetries are the same as in flat space: the only component different from zero can be taken to be 
antisymmetric in [fip]. Now though eq. (|49p gets an extra piece proportional to the Riemann tensor 
from the commutator of covariant derivatives: 

V.V^V^V^ = V [M (V^V p] V CT 0) + R e wa d e 0d u( f>. (50) 

The first term can be integrated by parts as before while from the second we see that it is possible 
to trade two external derivatives in the vertex for a factor of R. However, the minimal number of 
derivatives with n fields ip is still 2n and because of that a term of the form (di/)) 4 (d 2 i/j) n ~ 4 R cannot 
be generated. 
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B Calculation of the bispectra 



At leading order in slow roll the field 7Tfc(r) and its derivatives, which must be inserted in (|25p in 
order to compute the bispectrum, read in conformal time 

^-'^kr- it{T> ~ 4 ^w" r • ^-'^mkr^ - <5i) 

Hence, from all the modes in (|25|) there will be a factor 

' e i(k 1 +k 2 +k 3 )r r^'j 



which is common to all interactions. Additional factors for the different fields in each interaction 
must be accounted for according to the following rules: 

d 

vr — >l-ikiT, 7T — ► -Hkfr 2 , n — > H 2 kfr 2 (2 + ik iT ) , > -iHk iT . (53) 

a 

Since we are looking at interactions with six derivatives, there will also be an overall factor of H e . 
Recalling the definition A$ = ^i/ 2 /(4eMp 1 ), the generic bispectrum will thus look like 

5M H A 2 f° 

FMi(ki, k 2 , k 3 ) = -i 12g ^- 2 fc3fe3fe3 J &t P Ml {ki,T)e lktT + c.c. + permutations, (54) 

where the function P(ki,r) contains all the appropriate additional factors for the interaction as 
specified above, divided by £f 6 r 4 (because of the factor of r~ 4 in (|25p and the factor of H e we 
extracted) . 

It is convenient to shift to the adimensional integration variable y = kfT, and introduce a dummy 
parameter /x in the phase (e ty — > e l ^ v ) to be set to one at the end of the computation. By doing 
this one can substitute r — > —(i/ht)-£- in the function P, which can thus be pulled out from the 
integration: all is left to do is therefore the simple integral of a phase, e w = —i/n- After all 
these steps the final result is 



6eM^ kfk%k$k t 



+ permutations , (55) 



where the factor of 2 comes from the complex conjugate (since the result is real). In other words, 
the specific momentum dependence of each interaction is contained in the operator P that must be 
applied to before setting /u, = 1. The operator associated to a cubic interaction (but the argument 
can be generalized to any interaction) with n space gradients is given by the composition of 
i) 3 differential operators according to the following substitutions: 

, kid . fki\ 2 d 2 .. fh\ 2 f„ h d \ d 2 



kt dfi ' \ kt ) dfi 2 ' \ kt J V h dfi J dfi 2 



ii) n operators —(ki/kt)-^, one for each spatial derivative; 

iii) one operator [(l/fct)^] -4 to account for r -4 . 
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The operators generated by the interactions studied in this paper are 



k\k\k\ ( ^ _|_ d 



h dyu 



2 + 



k2_d_ 



2 + 



d \ d 2 



A: 2 \ 2 d 2 



1-^-1- 

k t da \h) da 2 



k t dfi J d/j? ' 

i_ f^Y— 

fct d/z V fcj / d^ 2 



and 



fc 2 A: 2 A ; 2 [(^-fc3)V 2 + ^A 



fct djii) ( 1 fet d/i) + 2 ( 2 + k t dfi) ( 2 + kt d/i) 



for Mi, M2 and M3 respectively. All of them must be applied to and evaluated at fj, = 1. 
The calculation of the different bispectra is now straightforward given the relation 



d n 1 



(-l) n n\ 



(57) 
(58) 



il 

d^ 2 

(59) 



(60) 



and gives the results quoted in Eqs. (|26|) . ()27|) and (|28p . 



C A template that vanishes in the squeezed limit 

For the analysis of CMB data the enfolded template, equation (j3Tj) . is good enough for constraining 
the non-Gaussianity induced by the operator proportional to M3 in (|20p . Since the enfolded template 
is a linear combination of the equilateral and orthogonal templates, its amplitude can be constrained 
using the recent limits on these shapes, as discussed in equation ([35D . 

However, the enfolded template fails to capture some important qualitative behaviour of the 
physical shape as we will now explain. There are competitive constraints on the non-Gaussianity 
from Large Scale Structure observations |32j. These arise from the fact that the non-Gaussian 
corrections to the halo bias Ab have a characteristic scale dependence and one expects this scale 
dependence to be sensitive to the squeezed limit of the non-Gaussian shape. Both the orthogonal 
and enfolded templates go to a constant in the squeezed limit, lim r2 _^o r 2 r 3^0-^ r 2> r 3) = const., and 
one could thus expect some scale dependent effect in the halo bias induced by this non-Gaussianity. 
However, physical shapes generated by the single-field inflation models studied in the literature go 
to zero in the squeezed limit \23\ |T7]. Indeed this can be explicitly checked in [17] for the model 
from which the orthogonal shape was inspired, and for the model under study in this paper. This 
means that one expects the scale dependence of A6 to be negligible for these models, even if they 
have a large cosine with templates which do induce a scale dependence. 

This motivates us to introduce a new factorizable template which apart from having a large cosine 
with the physical shape also behaves like the physical shape in the squeezed limiil 13 !. Let us warn 
the reader though that it is not clear whether the templates introduced in this Appendix are useful 
for computing other quantities which are interesting for LSS studies, as they might be sensitive to 
specific limits where the physical shape and the template differ. 



13 A similar approach has been taken in Appendix B of [T7]. Notice that in the published version of the 
paper, the proposed template contains some typos, that have been corrected in the most recent version. We 
thank L. Senatore for correspondence about this point. 
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The new template will contain also monomials that go like k~ 4 : this will give us additional 
freedom to impose that it behaves as k^ 1 in the squeezed limit, ki — > 0. One way to build the 
template is to take a linear combination of the new monomials for which the leading divergencies 
kf 4 and /c ; -3 cancel. We do not have enough freedom to cancel also the kf 2 divergence, but we can 
at least require that the divergence is independent of the direction we approach the limit. This will 
allow a cancellation with the "standard" monomials. The only combination with such properties is 
given by 

16 A; 2 1 1 

Fi(fci, k 2 , k 3 ) = 4 + - I2T4 - T2T4 + c y clic P erms - • ( 61 ) 

y rC\ rC^rC^ tjfv^fX/^ 13 2 3 

Now we can combine this with the usual monomials diverging as k~ 3 and k~ 2 

F 2 (h, h 2 , k 3 ) = pp - — - + c y clic p erms - > ( 62 ) 

1 2 1 2 3 2i 1 3 

^l.^M = 727272. (63) 

The residual kf 2 divergence of F\, which does not depend on the direction, can be cancelled by 
a proper addition of F 3 . This gives a new template, besides the standard equilateral one, with 
the correct squeezed limit. We have thus a one-parameter family of templates going as k^ 1 in the 
squeezed limit: 

F = 4/ NL A| [aF x + F 2 + 2(1 + a)F 3 ] , (64) 

where a is a free coefficient which can be fixed by requiring the cosine with the physical shapes to 
be maximum, and A is some normalization which can be fixed, for example, such that the template 
equals the local one in the equilateral limi10. The equilateral template corresponds to a = 0. In 
table [2] we show the values of a that maximize the cosine with the shapes of different physical models. 
In figure [3] we show the form of the template that approximates the shape of the non-Gaussianity 
generated by the operator proportional to M3. 
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